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THE COMPRESSIBLE FLOW PAST VARIOUS PLANE 
PROFILES REAR SONIC VELOCITY* 

By B. Gothert and K* H. Kavalki 


SUMMARY 


In an earlier report UM No. 1117 "by Gothert, entitled '^erechnung 
kompressibler ebener Stramvngen lei hohen Unterschall- 
Anblasegeschwindigkeiten" (Calculation of Compressible Plane Flows at 
High Subsonic' Airspeeds) , the single-source method was applied to the 
compressible flow around circles, ellipses, limes, and around an elongated 
body of revolution at different Mach numbers and the results compared as 
far as possible with the calculations by Lamia and Bosemann. Essentially,- 
it was found that with favorable source arrangement the single -source : 

method is in good 'agreement with the calculations of the same degree of 
approximation by Lamia and B username. Rear sonic velocity the number of 
steps must be increased considerably in order. to sufficiently approxi- 
mate the adiabatic curve. After exceeding a certain Mach number where 
local supersonic fields occur already, it was no longer possible, in 
spite of the substantially increased number of steps, to obtain a 
systematic solution because the calculation diverged. This result was 
interpreted to mean that above this point of divergence the symmetrical 
type of flaw ceases to exist and changes into the unsymmetrical type 
characterized by compressibility shocks- 


SURVEY OF THE APPROXIMATION METHOD USED 


In report UM No. 1117, it is explained how by moans of successive 
approximations the veioc-ity field, around plane profiles at high subsonic 
speeds can be computed. But instead of representing the solution in 
closed form as practiced in the conventional methods up to now, a 
numerical method was developed by which a large number of approximating 
steps can be computed, such that the labor for the individual step Is 
not increased, with progressive degree of approximation. As a consequence, 
the restriction was removed of being able to compute, in general, only 

^'Dle kampressibie Sbrdmung um verschiedene' ebere Profile in Nshe 
der Schallgeschwindigkeit. 11 Zonbrale fur wissenschaf tliches 
Berichtswesen der Luftfohrtforscliung des Generalluft- zeugme ister s 
(ZWB) Beilin- Aalershof , Untersuehungen und Mitteilurgen Nr. 1471, 

January 6, 1945* 
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two individual approximations owing to the increasing ..amount of paper- , 

work Involved, On the other hand, however, it was taken into tiie 

"bargain that even at increasing degree of .approximation the u mosh width" 

of the numerical method is, in general, maintained, so that a small 

consistent residual error remains, "but judged from the model problems 

worked out, this error remains, on the whole, far below the error 

introduced by an insufficient number of individual approximations. 

\ 

The afore-mentioned report proceeds from the concept that the 
compressible flow oan be regarded as a special incompressible flow, 
where the entire flow field outside of the profile in the stream is 
covered with elementary sources and sinks, the yield of which is 
dependent in a simple manner , upon the magnitude of the local Mach 
number. Thus, if the velocity field around a profile were known, the 
elementary sources themselves would be known also. The departure of 
the compressible from the incompressible velocity distribution could 
then be explained simply by the fact that the additional velocities 
produced by the sum of the elementary sources are computed. By the 
method pursued in the subsequently described model problems, the 
elementary sources were determined first from the velocity field of 
tho incompressible flow, that is, the velooity field of the Prandtl 
rule, and with this first approximation formula for the source inten- 
sities, the variations of the velocity field computed. With this 
improved velocity field, the source intensities are then computed again - 

and from it the velocities until with progressive degree of approximation 
the additional velocities due to the sources tend toward a fixed value 
that represents the desired velooity in the respective point of the ► 

compressible— flow field. 

✓ 

As the calculation of the elementary sources for given velocity 
is comparatively simple, it is important that the velocity field, which 
is produced by the sources distributed two dimensionally over the 
entire flow field, be computed in the simplest possible manner. For 
this purpose, the incompressible comparative flow covered with sources 
is conformably transformed on the circular flow, where the individual 
elementary sources do not vary in yield. The two dimensionally distri- 
buted sources were combined into concentric rings or, for even rougher 
approximation, into single-point sources which were disposed in suitable 
manner at such points where the maximum yield .of the sources is to be 
expected. Even this very rough method proved surprisingly satisfactory 
by suitably arranged single souroes for defining the velocities in the 
proximity of maximum profile thickness as verified by a comparison with 
the known calculations of Lamia and Busemann with the substantially more 
accurate source ring method. 

* 

The present report deals with a number of model calculations for 
profile flows worked out by the single-source method. The number of 
single sources in each flow gradient was fixed at twelve, the choice 
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of arrangement being such that the single sources disposed along a 
stream filament are of equal yield as much as possible, The .calculation 
vas further simplified by postulating symmetrical flows, Unsymmetrical 
supersonic flows with compressibility shoe Is are thus ruled out before- 
hand. 

One particular advantage of the single— source method cons jets in 
being able to see at once whether the chosen number of approximation 
steps already complies adequately with the continuity condition of the 
compressible flow or whether the method diverges and so indicates 
that, at the particular Mach number for the chosen profile, symmetrical 
flow free from compressibility shocks is no longer to be expected. 


II. VELOCITY DISTRIBUTION ABOUT A LINE (CRESCENT) (d/l = 0.10) 
1. Mach Numbers with Convergence of Computed Flows 


The chosen single-source scheme (system I, appendix) for a crescent 
(d/l = 0.10) at M«= 0.77 is represented in figure 1. The 12 single 
sources are arranged so as to lie particularly close together near the" 
point of maximum profile thickness in order, to be able to better account 
for the high yield of the sources expected here. But at greater distance 
from the profile, the sources are spaced farther apart because the total 
yield, owing to the already subsided increases of velocity and the 
reaction of these outer-lying sources to the velocities at the profile 
contour and on the areas of great source density, is substantially less. 
The additional velocities due to the single sources were computed in 
starting points disposed as much as possible at equal distance from the 
adjacent single sources in order to minimize the error due to the point 
concentration of the elementary sources. The represented source diagram 
varies with Increasing Mach number so that the individual sources move 
farther outward into the flow. This outward travel of the sources, 
desirable from the physical standpoint, offers mathematically an 
alleviation or facility by reason of the fact that after application of 
the Prandtl rule, the source locations in the plane of flow, conformally 
transformed in a circle, are coincident for all Mach numbers so that 
the additional velocities induced by the sources in the plane of the 
circle can be computed with the aid of the same Influence factor. 

The result of the approximate calculation with the previously 
ds scribed system of sources and starting points at = O .77 io 
illustrated in figure 1; the stream density (pv) is plotted against 
the M^ch number for different groups of starting points; the correct 
value of the stream density striven for by the approximate calculation 
is given by the adiabatio curve which at sonic velocity Mt * 1 has 
its maximum as sign for the greatest possible stream’ density. At point 
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Mt » M+oo, which corresponds to the state of the flow ■velocity, the 
adiabatic curve is tangent to the stream density curve of the Prandtl 
rule. On the Prandtl curve, all points of .the approximate calculation 
by the Prandtl rule are IcJCfSftpd; they are identified in the diagram as 
0^ h approximation, suoh as the starting point’ 11 at the maximum velocity, 
etc,, for example. The difference of the stream density in the particular 
approximation points relative to the accurate adiabatic curve indicates 
then the extent to which the stream density by Prandtl rule is still 
too great and to which the width of the stream filament passing through 
the respective starting point is still too small to fulfill the condition 
of continuity of the compressible flow. This difference represents, 
therefore, at the same time a measure for the additionally to be applied 
source yield, which is to force the particular stream filament apart 
to the correct width. On computing the' additive velocities by the sources 
which are dimensional according to the velocity field of the Prandtl 
rule, (0*& approximation) the velocity field of the next higher is 
obtained or, in this case, of the first approximation. It is true that 
the points of the first approximation found this way lie already much 
closer to the adiabatic curve than to the Prandtl curve; but it is seen 
at once that the second degree of approximation,, especially at the points 
of maximum velocities, is still not adequate to approximate the adiabatic 
curve closely enough. Since the still remaining difference relative to 
the adiabatic curve is immediately apparent from the diagram even for 
the points of the first approximation and the additive yield is thus 
known, the second approximation and also the higher approximations can 
he calculated at once. For the example in question, it is seen (fig. 1) 
that- about five approximative steps are necessary to approximate the 
adiabatic curve accurately enough starting from Prandtl ’s curvo. 

The number of steps required is, of course, less at small Mach 
numbers of ’ flow. Thus, figures 2 and 3 show, for example, that at 
Moo* 0.75 for the same profile only about three, and at tea = 0.70 
only about two steps are needed to attain the same degree of approxi- 
mation as at Moo» 0,77 with five steps. On the other hand, the 
number of steps required at higher M increases enormously. From all 
these approximate calculations reproduced in figures 1 to 3 j it is 
apparent that in the middle section (starting points 11 to 71) the 
Velocities in respeot to the Prandtl curve are still considerably 
increased; while at the starting points 12 to 32, farther upstream, 
this increase has already subsided appreciably. At the points 13 to 23, ■ 
farthest upstream, it has practically disappeared so that, on the whole, - • 
the velocity distribution along the contour decreases substantially with 
Increasing Mach number , ' 

2. Mach numbers with Divergence of the Computed Flows 

The approximate calculation at M» = 0,80 and M» =» 0.82 
was carried out by the same scheme and the results plotted in fig- 
ures 4- and 5. In spite of the much larger number of steps (up to 10), 
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no satisfactory approximation . of the adiabatic curve was obtainable. 
While at Moo - 0.80, the ninth and tenth step’s appear to approximate 
the theoretical curve to some small extent, divergence plainly 
noticeable at Moo = 0.82 in the fifth step even as is evident in 
the starting points 11 and 21. This divergence seems to be traceable 
to the starting points which extend into the range of supersonic speed 
(M+ > 1) . The significance of the divergence is discussed at the end 
of the report after the results of approximate calculations for sev- 
eral other profiles are available for the discussion. 


3. Discussion of the Computed Maximum Increases of Speed 
Compared with Busemann * s Data 


Figure 17 represents the velocities in starting point 11 at 
various, Mach numbers from which the increasing number of necessary 
pjteps with increasing M is apparent. When the velocity in starting 
point 11 nears the velocity of sound, this profile requires about 
five to six steps by the Prandtl rule to assure sufficiently close 
approximation. So, while the approximate calculations available so 
far for slender profiles were, because of the enormous paper ’work 
involved, carried out only to the first, in a few cases to the second 
step by the Prandtl rule, satisfactory compliance with the continuity 
condition at those high subsonic speeds is definibely not to be counted 
on. This holds' true so much more because by this approximation method 
the individual steps always yield a little less than by the single— source 
method, because their source yield was not applied in full strength 
corresponding to the momentary velocity field, but merely in an approxi- 
mation, the quality of which is equivalent to the degree of approxima- 
tion of the whole calculation. 

From the velocities obtained in the different starting points, 
the velocity. at the -point of maximum thickness was then computed by 
extrapolation with the aid of the tangent condition (dv = — dn, 

r = radius of curvature of stream. lines } , which at the same time 
represents the maximum, velocity occurring, at the crescent contour 
and plotted in figure 19 for the individual .approximations in com- 
parison with the values by Prandtl and Busemann. It was found that 
the first approximation of the single-source method was in good 
agreement with Busemann* s data. Thus, in this example, the arrange- 
ment of the single sources was obviously fortunate enough so that the 
errors of the rough approximation in the zone .of maximum increase of 
speed are nearly obviated. (However, it should be pointed out here 
that in other cases the agreement is not quite as good as for the 
crescent, for example, on the ellipse. This defect in the quality of 
the approximation method is, however, already, evidenced during the 



6 


MCA TM No. 1203 


calculation in the very nonuniform .distribution of the yield of the 
source along the individual streamlines so that a correction with 
better distribution of the single sources would have been possible 
at any time . ) 

After exceeding the limit of the velocity of sound, figure 19 again 
indicates that the number of necessary approximation steps increases 
considerably. A substantially greater number of steps beyond the fifth 
approximation were calculated and from these the values toward which 
the velocity tends at starting point 11 for an unlimited number of 
steps, determined by extrapolation or else, as in the immediate proximitj'- 
of the point of divergence, by a refined method. It resulted in a very 
accurately defined limiting Mach number at which the method ceased bo 
converge. The thus obtained point of divergence 13 e3 definitely in 
the supersonic range at a maximum local Mach number of M* = I.O96. 


4. Effect of Dissimilar Source Arrangements 


The arrangement of the sources must proceed from tne point o'f 
view of placing the sources sc that the areas of high yield of the 
elementary sources are covered as uniformly as possible with single 
sources. Particularly, the single sources lying on a stream filament 
should be distributed for most uniform yield. Since at the start of 
the approximation, in general, only sue incomplete estimation of the 
stream densities is possible, unfavorable systems of sources may be 
encountered. Thus, in order to gain an insight into tlie potential 
errors introduced, the crescent flow was computed again in the vicinity 
of the point of divergence with a modified source system, which in 
figure 16 i3 contrasted with the old system. In the new system III, 
the single sources are moved nearer to the area of maximum increase 
of speed and disposed more closely; owing to thi3 closer coverage near 
the maximum speed increases, the two outermost sources of system I were 
omitted so as not to increase the number of sources much above twelve. 

The results of the approximate calculations with both systems of sources 
for the starting points in the central section closest to the profile 
are shown in figure 16; once for Moo = 0.77 within the convergent zone, 
of the method, then a little above the diverging point at = 0.30. 

At M 00= O.77 both source systems manifest distinct convergence and 
at Moo = 0.80 divergence. To be sure, it appeal’s that the now 
source system II already ceases to converge at a somewhat -lover • 
limiting Mach number; but this displacement should be extremely little . 

In the case of convergence , the velocities themsc-lvee, resulting for 
the particular starting points at high degree of approximation, differ 
little from each other for both source systems. 

From this comparison, it can be concluded that, while dissimilar 
choice of source systems may result in slightly different numerical 
values for the starting point velocities, thin error remains considerably 
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below that caused "by an. insufficient number of approxi mating steps 
(as for example, breaking off the calculation after the first or second 
step beyond the Prandti rule). Even the fundamental characteristics 
of the computed flow, such as convergence and divergence of the solutions, 
remain completely unaffected by the type of special source system. 


Ill . TELOCITY DISTRIBUTION OYER ELLIPSE ABU ELONGATED 
PROFILE'' (d/7, = 0.1Q) 


As for the crescent, the compressible flow past an ellipse and a • ‘ 
spindle-3haped profile of the same d/l = 0.10 was computed and the 
results for several Mach numbers illustrated in figures 6 to 11. The 
source system I was used for the spindle-shaped profile while for the 
ellipse, a different arrangement of sources (system II) was chosen 
where the single sources extend, farther upstream or downstream so as ' 

to secure a better adaptation on the more complete velocity distribution 
of the ellipse. These profiles' also evince a point of divergence 
beyond which the applied calculation method ceases to converge.' This 
point lies at different Mach number of the basic flow, depending upon 
the particular profile form. But it is always observed when the 
velocity of sound is locally exceeded in the flow around the profile . 

This characteristic is especially reproduced for the ellipse in figure 17* 
where the velocities ere compared to several Mach numbers in the ' ' 
immediate proximity of maximum, profile thickness (starting point ll). 

The maximum increases of speed at the contour Av^/vo, of the 
ellipse were also computed for the several degrees of approximation 
and plotted in figure 19 along with the values of Busemann’s approximation. 
At the Mach number of flow corresponding to the divergence point, a 
local increase of speed of 2 6 .k percent results, equivalent to a ' ‘‘ 

maximum local Mach number of Mt may = 1.08, The comparison with 
Busemann’s method of approximation indicates that the first step of 
the single— source method yields a little mere than the corresponding 
Busemann step. Even though fundamentally the single-source method 
must yield a little more for the source intensities than Busemann’s 
method, on account of the absent series expansion, the existent dif- 
ference appears, nevertheless, to be largely attributable to the fact 
that the chosen system of sources was not well enough adapted to the 
flow around the ellipse as evidenced by the nonuniform distribution 
of the yields of the source. However, since the observed difference 
amounts, at the most, to 0.5 percent of the parallel flow, a correction 
of the calculation by means of a modified source system was omitted. 
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IT i : .VELOCITY DISTRIBUTION ABOUT A CIRCULAR CYLINDER IN 
• COMPARISON WITH LAMARS CALCULATIONS 

As a further example for calculation by the single-source method,, 
the compressible flow around a circular cylinder was chosen. The 
Prandtl rule certainly represents no good initial approximation for 
the flow around a circle, because the assumptions for the applicability 
of the Prandtl rule, such as regards sufficiently small increases of 
speed, for example, are not complied with. For this reason, the 
ordinary incompressible flow was used as the initial, approximation for 
the ,-circle so that the respective sources to be applied must cover the 
entire stream density difference between the incompressible and the 
compressible flow. The result of the calculation in the vicinity of 
the divergence point is reproduced in figures 12 to 15 . At Mm = 0.40, 
about four steps afford a sufficient approach to the exact adiabatic 
curve; at M OT = 0,42, about five steps are already necessary and at 
M w = 0.45, the calculated eight steps indicate that at this Mach number 
of flow, no convergence can be any longer counted upon; at M«> » 0.50, 
divergence definitely prevails. The velocities in the characteristic 
starting point 11 computed at different Mach numbers clearly show 
(fig. 18) how the number of necessary approximation steps from. 1 to 2 
at Mod = 0.25 up to the divergence point Increases considerably and 
ultimately after exceeding the divergence point, no approximation 
satisfying the adiabatic equation can be found no matter how many steps 
are used. 

The maximum increases of speed occurring at the circle circum- 
ference were computed from the data of the single-source calculation 
and plotted in figure 20 for different degrees of approximation and 
Mach numbers. Where the point of divergence lie3 at around M«,s 0.447 
with a maximum increase of Av^/v^, = 1.70, the equivalent to a maximum 
local Mach number Mh^^ = 1,295. The values of the first to third 
approximation by Lamia (reference 1) were also included. The comparison 
indicates good agreement between Lamia* s data and those obtained by the 
single-source method. 


V. SIGNIFICANCE OF THE POINT OF DIVERGENCE 

In all. of the afore-mentioned calculations, it Wa3 found that the 
single-source method, once a certain limiting Mach number had been 
exceeded, was unable to produce satisfactory approximate solutions 
because the calculation diverged. Typical for the divergence point was 
the fact that the highest local Mach number at the profile circumference 
was only slightly above the sonic velocity. To illustrate: on the 
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ellipse (d/2 = 0.10) the highest Mach muster related to the divergence 
point was = 1.082, for the crescent (d/2 = 0.10) = 1.096, 

and for the circle, t'T f max = 1.295* At this state of affairs/ the 
question arises whether perhaps by the rough approximation assumptions 
of the calculation only the maximum local Mach numbers are shifted into 
the supersonic range while by the exact calculation the sonic velocity 
just represents the critical limit at the exceeding of which the conver- 
gence stops. On top of that, it is extremely important to decide whether 
the observed divergence merely represents some defect of the mathematical 
method of calculation or whether it is perhaps the evidence of a specific 
physical progress. 

To analyze the last problem, the single— source method was modified 
so that the calculation was reduced to the solution of a 12— term 
quadratic equation system with 12 unknowns. (This investigation 
is to be published in a separate report.) After this change, it was 
possible, for example, to calculate flows past profiles with high local 
supersonic fields of M + rQa2 = 1.5 and more. This therefore proved that 
the single— source method" contains no fundamental difficulty to push 
ahead into the area of high local supersonic fields. It was also found 
that above a certain limiting Mach number, equivalent to the divergence 
point, a solution joining on the subsonic flow no longer exists. There- 
fore, the reason for the unsatisfactory approximate solution for the 
compressible flow around circular cylinders, ellipses, and crescent 
above the point of divergence is due to the fact that above the limiting 
Mach number, the symmetrical flow pattern postulated by the calculating 
formula no longer exists, but rather is replaced by an unsymmetri'cal 
flow associated with compressibility shocks. 

Concerning the maximum local Mach numbers related to the divergence 
point, it is appropriate to reexamine figures 17 2 nd l8 with the velocity 
variation in starting point 11. It is seen that the point of divergence 
is reached when the velocity curve related to this point 11 is exactly 
tangent to the adiabatic curve. But this contact point is always in 
the supersonic range as is to be inferred from the ascent of the velocity 
curve for point 11 and the adiabatic curve. The divergence point could 
lie then exactly at the sonic velocity only if the velocity curve of 
point 11 with horizontal tangent approached the adiat '.tic curve . 

The question is then whether the ascent of the velocity curve of 
point 11 can be varied as far as the horizontal entry in the adiabatic 
curve by refinement of the calculation formula. The refinement of the- 
formula can, for example, be visualized such that the previously 
employed 12 single sources are consistently split up until, in the 
limiting case of the exact calculation, infinitely many. small elementary 
sources are involved; but this splitting up modifies in no way the sum 
of the yields of tho source. Even though this splitting may slightly 
modify the magnitude of the additional velocities produced by the sources. 
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that is, shift the total velocity curve from the first approximation 
somewhat, it is extremely unlikely that the curve up to horizontal 
entry would he changed as a result. It is therefore necessary to 
reckon with the fact that the computed local supersonic fields do 
actually occur and can he symmetrically removed without compressibility 
' shock. By the worked-out model problems, the highest local supersonic 
Mach numbers still obtained in symmetrica! flow are so much greater 
as the Mach numbers of the airspeeds related to the divergence point 
are smaller. This is exemplified in figure 21 where the highest local 
Mach numbers of the divergence point are plotted against the airspeed 
Mach number for the circle, crescent, and ellipse. At Moo = 1 the 
curve of the. highest local Mach numbers has the value v = 1, 

hence at this speed the increases of velocity exactly disappear. It 
means that at sonic velocity M«> =1 a symmetrical flow is possible 
only for the infinitely thin, flat plate in parallel flow which is 
consistent with the variation of the flow density curve with its 
maximum at sonic velocity. Figuro 21 further shows the Mach numbers 
of flow at which local sonic velocity is exactly reached. A comparison 
with the Mach numbers of flow related to the divergence point indicates 
that both differ very little, the difference becomes less as becomes 

greater. At M«= 0,80 the difference amounts to a mere = 0,030. 

So, as far as the application of these data to the conditions in aero- 
technics is concerned, it is quite iramateritd. whether the flow abandons 
its symmetrical character at the sonic speed limit or at the Mach 
number of the divergence point and changes to unsymmetrical flow with 
compressibility shocks. This result holds for the present only for 
the discussed model flows around the circle, ellipse, and crescent. 
However, it is not very likely bhafc the conditions will be materially 
different on other profile forms. Since on the ellipse and on the 
crescent of equal d/2, for example, profile forms with very dissimilar 
curvature variation are already involved. 

From this point of view, the. question of whether the previously 
computed symmetrical flows with local supersonic flow fields are stable, 
that is, whether they return to the initial flow attitude after minor 
disturbances, itself loses. importance . In the calculation itself, only 
the conditions for potential flows, such as continuity, freedom from 
vortices, etc., were contained so that the problem of stability would 
have to he treated additionally end therefore it is not ascertainable .. 
for the time being. - 


VI. HESUISB- 


From the results of the calculations by the method of approxi- 
mation on compressible plane flows with the aid of single— source 
arrangements, the following findings were arrived at: 
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1. The maximum increase of speed on the circumference of a 
circular cylinder computed by the single— source method and compared 
with the data hy Lamia indicated that the calculated maximum speeds 
in all comparable degrees of approximation ( steps 1 to 3 beyond the 
incompressible flow) are in -very close agreement by both methods. 

In the case of the flow around an ellipse, a 0,5-percent 
difference in airspeed resulted in comparison with .Bus emann'e cal- 
culation ( approximation step 1 ), which cbyicusly is attributable 
to the not very favorable arrangement of the single source. But the 
magnitude of the 9rror involved did not seem to justify the additional 
paper work involved in a calculation with an improved source distribution. 


The first step for the flow around a crescent beyond the Frandtl 
rule could be compared with Busemaim J s calculations; here also the 
agreement was practically perfect . 

The conclusion, therefore, is that the single— source method is 
suitable for a quick, numerically correct insight into the highest 
increases of speed accompanying the compressible flow past profiles. 

2. As the single— source method affords a clear view of the 
quality of the attained approximation to the adiabatic curve, it 3s 
possible at all times to decide whether the attained degree of approxi- 
mation of the calculation is already sufficient or whether the number 
of steps needs to be increased. In this respect, the model problems 
worked out disclosed that at low airspeeds two to three steps are 
sufficient, but that in proximity of sonic velocity, six to eight steps 
are necessary. Thus in this range, the calculations by Rayleigh, Lamia, 
and Busemann would also have to be continued up to substantially higher 
degrees if sufficient approximation to the adiabatic curve in the sonic 
speed range is to be attained. 

3. After exceeding a certain Mach number of airspeed, it was no 
longer possible to secure a satisfactory symmetrical approximate 
solution with any number of steps because the calculation diverged. 

Since in other cases it succeeded in computing by the same method 
subsonic flows with large local supersonic zones (^j. = rv, 1,5) , the 
observed divergence cannot be attributed to the imperfect suitability 
of the employed calculating method. The appearance of the divergence 
was, therefore, interpreted as that the symmetrical subsonic flow can 
be continued only up to the Mach number related to the point of diver- 
gence and after that, the flow changes into the unsymmetricai type 
characterized by compressibility' shocks. 

k. The highest local velocity appealing at the point of diver- 
gence lies a little above sonic velocity, the -more so the smaller 
the Mach number of flow related to the point of divergence* 
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On the circle, for example, the maximum local Mach "number is 

M+ max = 1 ‘ 2 9 , at the crescent (d/l = 0 , 10 ) Mt * 1 . 10 , . and on the 

ellipse (d/l = 0 . 10 ) M+„_„ » 1 . 08 , max . 

5. The Mach .numbers of air flow related to the point of diver- 

gence are only a little higher than the Mach numbers of flow at which, 
locally at the circumference, the sonic velocity is exactly reached for 
the first time: the difference in. Mach number is, for example, 

AM-os* 0.030 for the ellipse (d/l ■» 0.10), AM*, = 0.030 for the 
crescent (d/z =: 0 . 10 ), and AM W s= O.O50 on the circle. Thus for 
actual flight practice, it is puito immaterial whether the flow already 
changes into the uneymmetrical type on reaching sonic speed or at the 
Mach number of the point of divergence . 

6. The flows with local supersonic fields computed by the single- 
source method satisfy the conditions for potential- flows such as 
continuity, nonturbulence, etc. However, it yet remains to be proved 
whether the calculated supersonic fields are stable,, that is, whether, 
after a small disturbance, they return in the initial state. 


Translated by J. Vance r 
National Advisory Committee 
.for Aeronautics 
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APPENDIX 

METHOD OF CALCULATION 

1- Choice of Starting Points and Source Systems 



Since the sources to be applied 
0 on the crescent and especially at the 


p 5 ! spindle-shaped tody are predominately 
placed within a comparatively narrow 
range near the profile center, a . 
system of sources was chosen at which 
the three -source lines -3* = ^Constant 
lie close together (system of points l). 
For the circular profile and especially 
for the ellipse on the other hand, the 
area of the additional sources is more 
extended toward the stagnation point. 
Owing to the steep increase of speed 
at the ellipse aft of the profile nose, 
the maximum, of the additional sources lies for this profile in the forward 
half of the profile even at high Mach numbers (system of points II) . For 
the study of the influence of the change of the source arrangement, a 
second system was included for the case of the crescent, where, instead 
of the 12 sources extending far into the outer space, 13 closely spaced 
sources were chosen (system III). 


Starting Point and Source System I 
(Polar Coordinates in the Plane of the Circle) 
[starting point (r^., ) , Source point (r^, ^*) , 

r Q = Badius of image cix-cle | 


Starting points 

(m-v) 

Source points 
(hi) 

Starting point and radii of circle of source 

r |W r o 

1.075 1.242 1.436 1.659 2-216 2.959 3.952 

Starting point angles 

fly = 90 
= 74 
= 6o 

pV (1 1) (2 1) (3 1) (4 1) (5 1 ) (6 1) (7 1 ) 
(1 2) (2 2) 2) 

(1 3 ) (2 3) 

Source -line angle 
= 78.75 
= 67-5 
= 52-5 

hi (1 1) (2 1 ) (1 1 ) 

(1 2) (2 2) (3 2) (4 2) (5 2) (6 2) (7 2) 
(1 3 ) (2 3 ) 
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Starting Point and Source System II 


Starting points 

(nv ) 

Source points 
(kl) 

Starting point and radii of circle of source 
1.075 1-242 1.436 1.659 2.061 2.972 3.952 

Starting-point angles 
d v = 90.0 
« 670 
■ 45-0 

y-V (1 1 ) (2 1 ) (3 1 ) (4 1 ) (5 1 } (6 1 ) (7 1) 
(1 2) (2 2) (3 2 ) 

(1 3) (2 3) , • 

Source-line angles 
= 71-75 
= 67-5 
= 56.25 
= 38,0 

kl (l 1 ) (2 l) (3 1 ) 

(4 2) (5 2 ) (6 2) (7 2 ) 

(1 2 ) (2 2 ) (3 2 ) 

(1 3) (2 3) 


Starting Point and Source System III 


Starting points 
(fiVf 

Source points 
(k\) 

Starting point and radii of circle of source 
1.075 1-242 1.436 1.659 2- 06 l 

Starting-point angles 


d v = 90.0 

|AV(1 1) (2 1) (3 1 ) (4 1) (5 1) 

= 78.75 

(l 2 ) (2 2 ) (l 2) 

= 670 

(1 3) (2 3) (3 3) (4 3) (5 3 ) 

Source -line angles 


= 84.35 

k\ (1 1 ) (2 1 ) (3 1 ) 

= 73-12 

(1 2 ) (2 2 ) (3 2 ) (4 2 ) (5 2 ) 

= 61 . 88 

• (1 3) (2 3 ) (3 3) (4 3) (5 3) 


2. Influence Factors in the Cyclic (circular) Plane 

, With fx* = P,e^* denoting the local vector of a source and 

x = re~^ the starting point, the velocity induced in the starting 
point is 

^Throughout this report x and v are substituted for the 
German script x and v which were used in the original German 
version. 
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The contour conditions are satisfied "by adding to each principal source 
(Hqu) the source (Spqu) reflected at the contour circle- If the 
source arrangement is in addition doubly symmetrical, that Is, sources 
of equal strength in ±d* sources and in ±(ji - $*) sources, four 
principal sources each can be combined. For 7 = (i5* - d) and 
7 = ~(d* + d), is to be assumed positive, for 7 = jt-(d*-iS) 

and 7 - - jt + (d* - d ) , negative, so that 



In order to better assess the influence of a continuous source 
distribution by a system of 12 single sources, the latter are visualized 
as being replaced by a source-line system with two-dimensional source 
strength distributed along rays d* = Constant. Combining then the 
influence of one radius at the source points by averaging and taking 
the source strength along the interval as constant, the influence of a 
step-like source distribution along the ray d* = Constant that sub- 
stitutes the line superposition is closely approximated although in 
the practical calculati on, only single sources are involved. Unless 
the influence quantity ^-fy$ plotted for fixed d and fixed v* 

against t varies excessively, the mean influence Is equal to the 
influence in the center, or within the. framework of the distribution 
by powers of 1-075 in the intermediate partial points. In the vicinity 
of the extremes of the influence curves at t = 1 , however, and at 
interval increases at the fourth source ring the average value of the 
interval must be computed more accurately. Indicating the average 
formation by overscoring and including the factor l/2rc in the influence 
factor, the final d component of the Induced velocity follows as 
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3 . Conformal Transformation and the Ensuing Velocity Distortion 

The conformal transformation = f (£ ) gives the flow around 

the profile in the ik* plane from, the flow around the image circle 
j p | = r Q . The complex potential x = cp + i^ is 



F 


OdO circle 



Circle 


The flow around the circle | £ j = r 0 as contour is with t, = re^ 
and X = t, + ~ 



To determine then the velocity of a specific profile in the ik*-plane> 
only the distortion needs to he computed- 

a. Crescent .— If 8 is the edge angle of the crescent and 


(cL/l) 


±k 


0 . 10 /p. 
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the thickness ratio of the 
of the crescent gives 


given comparative profile, 


the geometry 


tan 



The conformal transformation is attained ly 


z,.. * 
ilc 

ft , 

n 

^0 - n 

! r 0 x 

( 

z ik* "j 

— - + n 
■‘•’0 

i + 1 j 
v° / 

J 

z it* _ (1 - s 2n ) 

+ i l2S n sin nj3 J 


r o t , Q 2n 0 ,,n 

1 + S - 2fe cos Zip 

with 


i_ . , 

Se ip = r ° 


yZ- + 1 
- 1 o 


The quantities S and |3 
circular coordinates 


can he represented as functions of the 


ft 2 _ 


1 r ^ + 1-2-21 


cos d jL 1 


\ r 0 ^ 


1 + 2™ cos i3 zL 


tan 2 ( ,5/2) 


3 = — 


arc tan 2-E- sin $ 


JL*& 


• o 


arc tan 


2r/z 


(£) 


x/2 


fA 0 ) 


•2 


1 
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The profile parameter n is given "by the thickness ratio-. As 
exponent it indicates the multiplication of the angle it in uhe 
stagnation point of the circle- It is 


rcn = 2j t - 8 


or 


n - 2 


4 

- arc 

it 



The velocity distortion follows from 


as 



(4) 


6 


= 2 arc tan 


Vs - 1 



- 2 


arc tan 



± 1 
- 1 



The formula of the conformal transformation was so chosen that parallel 
flow prevails at infinity. Actually it Is 



= 1, 3 = O' 


. 2 — 

? 3in * 

sin 2 n^ 


-> 1 


l/2 


which conditions that 


n. 


20 


NACA TM No* 1203 


b^Elong ^ed body of revolution. -The conformal trans. 


is supplied by the formula satisfying 


dZ., * ' 

ik 


formation 


^ = L + /i , 1- K /r 0 \ 3 


r K 
:° r/r c 


JOS -5 + y_.J? (y) cos 3-3 


K > . 

" 7 — ) si 
r/r oS 


Bin d - ~y~ (y) sin 3d 


Since 


V2 2 + k 


parameter K is 


p 2 (i - E) . 

“'3 * 18 obtained, the profile 


2 1 - 


.1/ik 


Z/ik 


Ay 

1 - E __ V/Lk 

3 % + l r 


The velocity distortion folio,,, from the ik* - velocity at 


7 X - iv y 


ik ■ 1 


(?)] [ x + (i - * cry] 
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With 


N = 




- K) cos 26 + (1 - K) 2 


we get 



hence , with ( 3 ) 


w 2 = 


1 + (1 ' K)2 + 2 


T -~) (l - K) cos 23 


1 + 


2 (~)“ cos 2-3 


= A + B cos 2-3 + C cos 43 


where 
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c. Ellipse. 
transformation: 


A corresponding formula supplies the conformal 


with 



as profile parameter. 


The velocity distortion follows from 



hence, with ( 3 ) the velocity hecomes 


cos 


i ♦(?} - <JX cos “ 


2-a 4 = */2 x 1 + (^o/r) 2 L li^y^r-7^, 

1 ♦ (a/r) 2 1 + < 4 / l >ik 77-^2 


r = r r 


2 sin $ 


1 + (a/r 0 ) *” “ 2 ( a / r o) cos 


r; , ' ^ 

1 V.1 ik r”p /<•]>. *2 2, 

L | gin ’9 + \^JJ COB ^ 


4 . Determination of Source Strength for Assessing the Difference 

Density and the Employed Approximation 


Between the Exact Stream 


If 0 _ ~ k .~ is the stream density of the < 

pX a x 

M * _ m = “ (a = local sonic velocity, a- 

a * ; ™ a oo 

result of the adiahatic phase change 


compressible f 3 .ow, we get with 
F = critical velocity of sound) as a 


^ 

0oo ~ ML 


K “ 1 w 2( M X 


* 1 " “ M <* 


'Jc - 1 


with 


B B a 1.405. Considering only the linear terms in the interference velocxby 


)onformably to Prandtl and putting 
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2k 


gives with 


Mfc # _ ( £%* 

Mco* " Moo* 

1 - vr 2 _ 1 

•1. Mco "*"q 


9, 

with 

yielding 


‘.r ■ rn 


AM k * 

“» i 1 + — 


,N 


MccV V 


M» e 


2 a m ) „ ! + 1 eas: = A ^? 
vv 




« 1 + -^ “T~ = Ay^- + B 




M * 

CO 


A » v. B - v 

M* 

0 Pr ' 1 /Mk* ^ 

- 1 + -si 1 , 


9 


H 


■w 


(8) 


as lineai'ized Prandtl -stream density. 

Differentiation of (7) gives as slope - of the stream density curve 
at point Mq* " 

&0L 

k - 


_£c^_ 

dM k * 

"M? 


1. 


. M 2\ 0 k MoA 

7 ^ 1 8„ Mv* ^~=TiJ y ^ 


The Prandtl -stream density 0p r , therefore, represents the tangent to 

the stream density curve in the approach flow point, while the incom- 
pressible stream density at speed M k * 

‘ '■!. 

" M w * " X *' Mi ' ‘ {9) 

can be regarded as secant through the points M k * = 0 and M^* = Mco*> 

Fow.dn order to be able to secure an improved solution for the 
compressible flow with the aid of the approximated stream density 
% = e p r or 0 ik , respectively, the compressible flow in a flow tube 
is visualized as being replaced by a flow of approximated stream 
density fljj with identical velocity field, in which the additional 
displacement is produced by an additional source distribution* 

(Compare UM 1117, p* 3)* If doy. is the width of the flow tube at 
a given point, ** 
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ffik _ fate e _k 

r 0 Vco ' r o Q oo 


will be valid for the flow along the flow tube- The approximation 
flow with additional sources is -to show the same velocities as the 
compressible flow* Due to the source distribution, the flow then 
differs from point to point; 



i’oVflo 


d R 
j^O ®co 


Both 0^ and 9-§ are functions of the velocity At 

infinity, the flow is the same in both oases, due to 0jj = 9 «,* The 

sum of all volume sources to be placed upstream then is (compare 
DM 1117, p* 28) 



If one introduces the constant flow dy._, which is calculated, say,, 
in the center section of the doubly -symmetrical profile, 




/ o.n- K 6 


s^O 


0o 


with Cw*{ * — ) = 

center section 


( 10 ) 


is valid. Cj\j* denotes, therefore, the relative error of the approxi- 
mated stream density in contrast to the exact stream density for the 
compressible velocity M*. 

a. Incompressible basic flow (flow around circular profile ) * - The 
curve of the source strength to be disposed about the circular profile 
in the incompressible flow is, therefore 


dE ilc = C-rv* t “ 


"o 9c y 


ik 




center section 


with C-y c * = 


'ik 


- 0 


k 


do 1 ) 


The width of the flow tubes at the circle in the center section is 
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b. Prandtl basic flow (crescent, spindle-shaped body, ellipse) . - 
For their profiles the Prandtl approximation is considered as comparative 
flow and the profile analyzed in the affinely distorted ik*-plane 
instead of in the compressible plane. With. denoting the velocity 

at the profile of the ik*-pl.ane, as_it is attained 'from the cyclic 
flow free from sources by conformal transformation, the velocity of 
the Prandtl approximation 2^*^ on the profile of the k- piano 
obtained by the transformation 


„ i + ^ ^ 

Ho* v m “* ■ y 


(id 


corresponding to the linearization of, the stream density in (8). The 
related stream density, therefore, is 

e ?ro ^ + 1_ / Hk*o %k* o 

0co~ ~ d 2 \mJ J ~ m m * 

e ik* M ik* 

with ~r — ~rf— denoting the incompressible stream density related 

& CO ^*00 .. , " 

to 2^*0 in the ik*-plane, the transformation of the stream 
density follows as 


~ ®ik* 


(12) 


low, if the Prandtl flow M * in the compressible plane corresponds 
to the flow M ik * 0 in the ik*-plane, the Prandtl i'low 
provided with additional sources replacing the compressible flow 
corresponds to a flow in the ik*-plane fitted with sources, the 
Velocity of which is equal to the velocity transformed by (ll) 


M iH* _ 
'Mco* = 


1 /%* 

*•**«'• V 

\M <# 



The throughflow of the Mq k * - flow with sources at a point of flow, 
tube width dn ik *, is then, since by (12) 0j_k* " Prandtl is the 

related stream density, 

^Ik* ®Pr 
r 0 9 a, 


%k 
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Oiling to the source distribution this throughfiow varies from cue 


do, v * 


* = 


place to another. An infinite distance d\jf ik 
the sum of the source strengths in the ik*~plsx&e 


IK 


•o. . 


hence , 


~ dn, ;k * /§p r 


fej£^\ 
“i kV 


For their profiles, however. 


Hk' 


on?- 

n co 


czyv* dn k 

compressihle flow dnj^^ = ah k 3 W by reason of the continuity. 
Therefore, the source strength to be applied in the ik*- j.laafc is • 


and in the corresponding 


^ik* = 


to ik* 


Pr 


- Q 


- 1 " Cp. r» 


^ik* a k 


r - ** I , with C-qv> 

y cq / * - 


* a 


9 P x- ~ e k 

0lr 


:io") 


while. by (10) the k-plcne follows at 


'dn k 0v 
*2k* = Cp r * ' - - K 


£ 1 = h dE ik * 


The source strength per width of flow tube remains, therefore, unchanged 
by the transformation. On comparing (lO n ) with (10 1 ) it is seen that, 
owing to (12), the same source strength dF lk * had to be placed in 
the ik*-plane, if a compressible flow was to be computed by moans, 
of additive sources, which is to flow through" the flow tube of 
the ik*-plane and there attain the velocity M.., * of the incom- 

° ±K 

pres si ole flow with sources- 


5- Iteration Process 

Since in the formula for computing the source strength the 
unknown compressible velocity Mv* enters on t2ae righthand side of (10), 
the source strength must be calculated step by step. With v 

as the velocity in the starting point pV of the comparative flow 
used as basis, that is, of the incompressible circular flow or Prandtl 
flow around one of the thin profiles the relative error of the related 
approximation stream density 0jj(M k * o ). in respect to the actual , . 

stream density 9 (M k *~ 0 ) is computed at the points i^r* 0 \j.V ul ’ fcer 

which (10) affords a first approximation for the looked for source 
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sum 


^ r Q V co 


= Generally,, if ( y 

' “1 


are the velocities 


of the (n^l)th approximation, ; thee source strength sum of nth approxi- 
mation is 


%J n> 


pV 


SP 


>“1) 




(n-l) 


0 k M},* 


(n-l) 


1 

- 1 


(13) 


with 


ap 


(n-l) 


Ap 


9 m 


J.V 


i i 


where ®k\Mk* / is preferably t alien from a diagram, and Ap 

represents the width of the stream tube in midsection» It is 
Ar __ 


Ap 


and 


a _ 5 respectively. Since the outermost source 
o v “ y " 2 

ring must include all the sour'ces existing the outer space, we get, 

for example, for the residual, source ring y- = 2*952, if ~ = 3-298 

' o v. > ° 

Ar > 


represents the limiting ring of the influence ring j related 

VO/fi 


to 


6 


A 7 = 


yik* 


(27) 




(K7) 


# 


+ 0.38 ~~~ (r 7 ) 


1-38 (r-^) 


Eik 7 


(E7) 


where O.38 indicates a reduction factor for the influence of the 
outer space. The individual source strengths Ejj^W follow then 

from the difference of two source sums in adjacent points on the 
same ring, since they relate to the same flow tubes. The additive 

velocities of nth approximation induced by the sources 

in the starting points pV are obtained by means of the influence 
factors ■ fg* at the circle, by neglecting the r-component and talcing 
the -q -component equal to the total velocity. (Compare DM Ul-T, p- 3-) 

fFor thin profiles for which the velocity differs little from the 
approved flow, the exact stream density curve (7) can, in vicinity of 
sonic velocity, be closely approximated by a parabola, (compare 
Oswatitsch) . On expanding the bracket in power of AM^* = - 1 and 

discounting all but the squared terms, we get d\ z 2 1 - ~g-~- (1 - Mk*) • 
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Tile influence factors for the Telocity 
the circular flow are, because of (2) 



simply 


of the source kX of 


uV 

~ 1 '9 


# 


The total effect of all sources in the starting point 11 V is obtained 
then by summation oyer all sources kX 




jj.v ~ (n) 

M «' ;; r KX \x 


(14) 


For the profile flow the influence factors (2) must be, in correspondence 
with the conformal transformation, divided by the distortion quantities 

(4, 5, 6): 


# 


_ j3__ 


kX 


F< iv 


and the additive velocities according to Prandtl distorted corresponding 
to (ll) 


l/at* 

l & 


(n) 


JdV 


= v-\* 


x — 

k! 


-ft'V Z (n) 
"kX kX 


(14) 


with this velocity connection the total velocities of the nth approxi- 
mation are: 


Mfe *(n) = _ (15) 

If to the initial velocity there corresponds the stream density 

of the basic comparative flow = 9 ik(%-*^) or e p r (^* 0 ) * 

respectively, then to the improved velocity " there pertains 

already a stream density nearer to the correct 0 -curve \ is 

calculated from ( 13 ) by solution of tlfe equation with respect to 0^., 

and by substituting for all magnitudes on the right side the values 
that correspond to the new approximation 
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d k 


(a) . 




n; 


1 + 


1 A 1+1 A. 


(n) 


ap 


(a) 


with 


a(i 


(a) 


Connecting the points Mv*^ , <9, 

^ 9 i£ * 04.^ gj.v©a 

8 Se ? ? which at convergence of the method approach the 

thra£eady y attained°?f S ° f approximation dre a picture of 

already aptained approximation of the 0 -curve. 


% 


(a) 


Mk* 


■(a) - 


a 


'co 




(16) 


for each starting point gives 
















Figure 2.- Line (crescent) (d/j = 0.10), M m = 0.70. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system I) . 
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Figure 4.- Line (crescent) (d/j =0.10), M TO =0.80. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system I). 
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Figure 5.- Line (crescent) (d/l = 0.10), M*, = 0.82. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system I) . 













Figure 6.- Ellipse (d/i = 0.10), M ro = 0.80. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system E.) 
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Figure 7.- Ellipse (d/Z- = 0.10), = 0.82. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system II). 


NACA TM No- 1203 














o.&s 


Stream density 
p v I 1 




— i — 
r Starting 

V , nrvi n+ 1 9. 

V 

i 

c 

Pi 

randl 

.l’s I 

•ule- 


>*? 
> . 

\ ' 
S 

\1 _ 


T 





_ 3 2 


x 

V. 

V 

N 

5k- 





— 


-Adi 

abat 

n 

LC CU 

rve 

> 







- Pa 

ralle 

lflo 

w M c 

* 

& 










' M 



1 

ach numt 

i 

>er a 


,x v K 
7 C -s 


Stream density 
-/.v 1 i L _ 




7~ 

j . 

Starting 


« ‘ 

Prand 

tPs 3 

-ule- 



/ 

% 

X 

pom 

i 

t 13 

T 




* 

23 

'xvrw. 


- -v 

**P 

/ 

it 

.3* 

n. 







-Adi 

abat 

n 

ic ci 








irve- 


-Paj 

•alia 

flo\ 











-* 

- Me 


tch n 
— 

umb 



er M 


x v K 


086 0.66 0.90 0.92 0.9* 0.96 0.9* 100 102 10* 1.00 0.6 6 0.66 490 052 0.9* 0.56 0.56 1.00 1.0Z 1.0* f.06 


Figure8.- Ellipse(d/J =0.10), M m = 0.838. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system n). 
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Figure 9.- ■ Spindle-shaped body(d/l = 0.10), M«>=0.76^ 

Velocity and stream density for various approximations 
to the adiabatic curve (source system I). 










Figure 10.- Spindle-shaped body-(d/J = 0.10), = 0.78. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system I). 
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Figure 11 



Spindle-shaped body (d/j = 0.10) , M 00 = 0.80. 

Velocity and streana density for various approximations 
to the adiabatic curve (source system I). 
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Figure 12.- Circle. M TO = 0.40. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system IT). 
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Figure 13.- Circle. Moo= 0.4-2. 

Velocity and stream density for various approximations 
to the adiabatic curve (source system n). 
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V elocity and stream density for various approximations 
to the adiabatic curve (source system H). 
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Figure 15.- Circle. Moo =0.50. 

Velocity and stream density for various approximations 
to the adiabatic curve (source systeml). 
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Figure 16. - [Title almost wholly illegible.] 

Various source systems 

on the line (crescent) (d /l = 0.10). 
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Figure 17.- 

Velocity at starting point 11 for various degre' 
approximation for the line (crescent) (d /l = 0 
and the ellipse (d/i = 0.10). 

es of 
.10) 














-£~ 

CO 



0.4 OS 0.6 0.7 OS 0.9 10 U 1.2 1.3 1.4 IS 1.6 

Figure 18.- Velocity at the starting point 11 for various degrees of 
approximation for the circular cylinder. 


i t r I r 1 


me A TM No. 1203 




Figure 19.- Maximum excessive velocity * -~-1 on the contour according to various 

V 


approximation calculations for line (crescent) (d/l = 0.10) and ellipse (d/l =0.10). 
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Figure 20.- Maximum excessive velocity on the circle contour 
according to various approximation calculations. 













